The Origin of Anomalous Low- Temperature Downturns in the Thermal Conductivity 

of Cuprates 
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We show that the anomalous decrease in the thermal conductivity of cuprates below 300 mK, 
as has been observed recently in several cuprate materials including Pr2-a;Ce2,Cu07_5 in the field- 
induced normal state, is due to the thermal decoupling of phonons and electrons in the sample. Upon 
lowering the temperature, the phonon-electron heat transfer rate decreases and, as a result, a heat 
current bottleneck develops between the phonons, which can in some cases be primarily responsible 
for heating the sample, and the electrons. The contribution that the electrons make to the total 
low-T heat current is thus limited by the phonon-electron heat transfer rate, and falls rapidly with 
decreasing temperature, resulting in the apparent low-T downturn of the thermal conductivity. We 
obtain the temperature and magnetic field dependence of the low-T thermal conductivity in the 
presence of phonon-electron thermal decoupling and find good agreement with the data in both the 
normal and superconducting states. 



PACS numbers: PACS numbers: 72.14.eb,74.72.-h 



I. INTRODUCTION 

In recent measurements made of the thermal conduc- 
tivity K of high Tc cuprates at milliKelvin temperatures, 
an unexplained decrease in the electronic component of 
thermal conductivity with decreasing temperature well 
below 1 K has been observed (we will refer to such de- 
creases as "downturns" ) . The most dramatic observation 
of a downturn in k/T occurred in the field induced nor- 
mal state of optimally doped Pr2-xCexCuOT-s (hence- 
forth PCCO), where the electronic contribution Kei/T 
was found to decrease below 300 mK and vanish at the 
lowest measured temperaturesi. This was in startling 
contrast to electrical conductivity measurements made 
on the same material (and to NMR studies^), which 
showed apparently normal metallic behaviour over the 
same temperature range. These measurements appeared 
to show an unexpected violation of the Wiedemann- Franz 
law at very low temperatures^ , with heat conductivity be- 
coming much less than charge conductivity. (Note that 
the main effect, observed above 300 mK, is a violation 
whereby heat conductivity is greater than charge con- 
ductivity by a factor of 1.7-1.8.)^ 

Downturns have been observed in other cuprate ma- 
terials in both the normal and superconducting states. 
For example, data on non-superconducting, overdoped 
La2-a:Sr2;Cu04 and on optimally and overdoped samples 
of superconducting La2-2:Sr2:Cu04 and YBa2Cu307_5 
(henceforth LSCO and YBCO, respectively) have showed 
similar downturns'^"'^^. In the d-wave superconducting 
state the downturn is associated with a decrease in the 
nodal quasiparticle contribution to k/T at a temperature 
well below that corresponding to the impurity scattering 
rate. This cannot be understood within standard dirty 
d-wave transport theorji. 



In this article we show how the observed downturns can 
be attributed to a decrease in the rate of energy trans- 
ferred from the phonons that provide the heat current to 
the electrons being studied (from now on we will use the 
word "electron" to refer to either a Landau quasiparticle 
in the Fermi liquid or a nodal Bogolubov quasiparticle 
in the d-wave superconductor). As a result of the re- 
duced heat transfer rate, the phonons no longer come 
into thermal equilibrium with the electrons at very low 
temperatures and the electronic heat current becomes de- 
pendent on the phonon-electron heat transfer rate, and 
thus decreases rapidly with temperature. With the origin 
of the downturn thus understood, the low-T thermal con- 
ductivity data for cuprates in the superconducting state 
are consistent with the predictions of dirty-d-wave the- 
ory. Similarly, the thermal conductivity of normal state 
PCCO at very low temperatures exhibits a constant lin- 
ear term in k/T versus T down to the lowest tempera- 
tures, consistent with the behaviour measured above 300 
mK, with no downturn. A detailed discussion of the ex- 
perimental results in a variety of materials will appear in 
Ref. i. 

We develop a theoretical model of the temperature 
and field dependence of the downturn by calculating the 
phonon-electron heat transfer rate and using it to deter- 
mine the measured heat current and the thermal conduc- 
tivity. In summary, our results (i) give the correct tem- 
perature dependence for the observed thermal conduc- 
tivity (ii) give the correct order of magnitude (which de- 
pends on the electron-phonon matrix element) of the ob- 
served effect (iii) give the correct dependence of the ther- 
mal conductivity on magnetic field in the vortex state. 

In Section II of this article we describe the role that 
phonon-electron heat transfer plays in low-temperature 
thermal conductivity measurements and explain the ori- 



2 



gin of the downturn. In Section III we show the calcu- 
lation of the relevant electron-phonon heat transfer rate 
in the normal state and in the vortex state of a d-wave 
superconductor as a function of temperature and mag- 
netic field. In Section IV wc will compare our results 
with data. In Section V we discuss some implications of 
our results. 



II. THE EFFECT OF PHONON-ELECTRON 
HEAT TRANSFER ON THE OBSERVED 
THERMAL CONDUCTIVITY AT LOW 
TEMPERATURE 

In order to study the effect of phonon-electron heat 
transfer rate on the low-T thermal conductivity data, we 
use a simplified model (Fig. ^) of the experimental con- 
figuration. A thermal current, Q, is carried into the sam- 
ple by phonons and electrons. The "contact" resistors, 
labelled Rei(c) and Rph(c)i correspond to the thermal re- 
sistance encountered by the electronic and phononic heat 
currents before entering the sample. Within the cuprate 
material, the thermal resistance to the electron heat cur- 
rent, i?ei, and to the phonon heat current Rph are as- 
sumed to be described by the standard theory, i.e. by 
Fermi liquid transport theory for the normal state and 
dirty-d-wave transport theory for the superconducting 
state. The temperature and field dependent phonon- 
electron thermal resistance Rei-ph, describes the heat 
transfer between phonons and electrons in the sample. 

The model of Fig. ^is clearly a crude approximation 
of the actual experimental configuration, but it turns out 
to be sufficient to describe the effect of phonon-electron 
heat transfer on the low-T thermal conductivity mea- 
surements. We will here briefly mention some of the 
differences between the model shown and the actual ex- 
periment. According to Fig. ^ the temperature differ- 
ence between phonons and electrons exists only at the 
hot end of the sample and the transfer of heat between 
phonons and electrons occurs before either the phonon 
or electron heat currents encounter any resistance from 
the sample. In reality, a temperature difference between 
phonons and electrons can be present everywhere in the 
sample and thus heat transfer can occur throughout the 
sample volume. The temperature difference that drives 
the heat current through R^i-ph actually depends on po- 
sition within the sample. However, since all our results 
below are calculated to linear order in the temperature 
difference between electrons and phonons, we can deter- 
mine the total phonon-electron heat transfer rate below 
as if a constant temperature difference existed through- 
out the sample volume. Also, in Fig. ^ we have shown 
only one set of current contacts, at the hot end of the 
sample, although a similar set of contacts must be present 
at the cold end of sample in the experimental configura- 
tion. One could slightly improve the model of Fig. 
by including a second, equivalent, set of resistors: Rei{c)y 
Rph(c) and Rei~ph at the cold end of the sample (the right 
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FIG. 1: A simple picture of the experimental configuration 
for thermal conductivity measurements. The thermal resis- 
tance to phonon and electron heat currents that occurs before 
the current enters the sample are represented by Rph{c) and 
Relic), respectively. The thermal resistance to phonon and 
electron heat flow through the sample are represented by Rph 
and Rei, respectively. The electron-phonon heat transfer rate 
is associated with R^i-ph- Temperatures at different positions 
have been labelled. 



side of the figure). This would account for effects of the 
thermal resistance encountered by the phonon and elec- 
tron heat currents after leaving the sample and would al- 
low for a temperature difference and resulting heat trans- 
fer between electrons and phonons to occur at the cold 
end of the sample as well as the hot end. We have deter- 
mined the thermal conductivity for this modified model 
and have found that the effect of the modification to our 
main result below (Eq. O is the same as that of simply 
doubling both the value of i?ei(c) and the effective sample 
length over which phonon-electron heat transfer occurs. 
Since we can only make rough comparisons of Rei(c) (as 
we discuss in Section IV) and the sample length (see foot- 
note, Ref. 0) with measured quantities, this modification 
does not significantly change any of our results below. We 
include these factors of two in comparisons with data be- 
low but will avoid needlessly cluttering the figure shown 
by including additional resistors. 

As indicated in the introduction, the electron-phonon 
thermal resistance Rei-ph plays a central role in the de- 
termination of the low-T behaviour of the apparent elec- 
tronic contribution to the thermal conductivity. In Sec- 
tion III, we calculate the electron-phonon thermal resis- 
tance for both the normal metallic state and the d-wave 
superconducting state and find that the result in each 
case is of the following form: 



(1) 



where fl is the sample volume, if is a constant, which 
depends on the electron-phonon matrix element, that we 
determine and n is between 4 and 5. For the remainder 
of this section, we will use Eq. ^ along with Fig. ^ in 
order to determine the low-T behaviour of the apparent 
electronic thermal conductivity. This is done to provide 
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a general explanation for the observed low-T downturns 
in both the normal state and the d-wave superconducting 
state before getting into the details of the calculation of 
Rei-ph for these separate cases. 

According to Eq. ^ Rei-ph varies more rapidly with 
temperature than any of the other resistors in Fig. ^ At 
sufficiently high temperatures, Rei-ph is negligible and 
the phonons and electrons at the hot end of the sample 
are in thermal equilibrium, so — T4. The thermal 
conductivity n can then be obtained from the thermal 
conductance G, where G = Q/(Ti - T2) = R~i^ + R~h- 
The result for either the Fermi liquid or the dirty-d-wave 
superconductor using the scattering rates dominant at 
low-T then has the usual form 



k/T ■ 



pj,2 (t»Td) 



(2) 



where the aT term is associated with the electron heat 
current and the I3T^ term is associated with the phonon 
heat current. We have introduced the characteristic 
downturn temperature Tp, which corresponds to the 
temperature at which R^i-ph = Rei, that is the relevant 
temperature scale for the downturn in the electronic con- 
tribution to K. 

At very low temperatures, Rei-ph becomes large, and 
the phonons and electrons in the sample become ther- 
mally decoupled. There arises a temperature difference 
between the electrons and the phonons at the hot end 
of the sample, T3 ^ T4, that depends on the relative 
magnitude of R^i^c) ^^id Rph{c)- There is thus some am- 
biguity in the definition of the thermal conductivity. A 
thermometer at the hot end of the sample could measure 
T3, T4, or some function of the two, say T' = f(T^,T4) 
depending on its relative sensitivity to phonons and to 
electrons. The thermal conductivity at low-temperature 
would be determined from a conductance given by G = 
Q/{T' — T2). One could consider many possible cases^° 
resulting from different forms of f{T^,T4) over the rel- 
evant temperature range. Here we focus on one limit- 
ing case, f{Tz^Ti) = r4 corresponding to a thermometer 
sensitive only to phonon properties. In this limiting case, 
an exact analysis of the configuration of Fig. Ogives the 
conductance 
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where 



A 



Rel- 



et— ph 



Rel-ph + Rel(c) + Rph(c) 



In this article we do not try to understand the nature 
of the contact resistances i?e;(c) and Rph{c) by looking 
into the details of the experiments, but instead deter- 
mine what relative values these quantities must have in 
order for us to be able to account for the observed low- 
temperature thermal conductivity downturns. To this 



end, note that when Rei(c) 
the thermal conductance 

Q 



becomes infinite, one finds 



T4-T2 



Rp 



R 



■el—ph 



1 

Rph 



(4) 



where again we have assumed that the hot end of the 
sample measures the phonon temperature T4. At T >> 
Td where Rei-ph « Rei, Eq- EJgives Eq. |21above. How- 
ever, at temperatures T << Tjj, where Rei-ph describes 
a high resistance to heat flow between phonons and elec- 
trons, the effective electronic contribution to the heat 
conductivity (the first term on the right side of Eq. ^ is 
much reduced. This is our picture of the low-temperature 
downturns. Alternatively, if the value of Rph{c) is taken 
to be infinite, the determined thermal conductivity shows 
no low-T downturn (this situation is discussed briefly in 
Ref. Hol) . To retain the simplest model of the exper- 
imental configuration that gives the correct low-T be- 
haviour, we can take R^i^c) >> Rph(c) over the entire 
temperature-range considered^i. 

Following this logic, we neglect terms of order 
Rph{c) / Rei(c) compared to unity in Eq. |21 We then ob- 
tain, after some rearrangement and the use of Eqs. ^^J^d 
|5]the thermal conductivity given by 



k/T = a- 



I r 

^ l+r(T/TD)"-i 



+ /3T2 



(5) 



where Td = (a/XZ2)i/(n-i) ^ is the length of the sam- 
ple along the direction of the current and r = i?e/(c)/-Re/ 
will be assumed to be independent of temperature. This 
is a generalization of Eq. ^ valid for finite i?e;(c)- At 
temperatures much larger than To the phonons and elec- 
trons in the sample are in local thermal equilibrium and 
Eq. [SI reduces to Eq. |21 At temperatures well below Td 
the electronic contribution is limited by phonon-electron 
heat transfer, so that k cx R^iLph due to the presence of 
the heat current bottleneck. Thus the electronic contri- 
bution decreases to its minimum value, kq — -l-r), 
as K — kq cx {T/Tjj)"' at low temperature. 

Eq. |5l along with its field-dependent generalization 
(Eq. I15|l that applies in the vortex state, is the main 
result of this article. If we use the values of n and Td 
that we calculate in Section HI, then Eq. Ogives a good 
description of the low-T data for downturns observed in 
the normal and superconducting states of cuprates. 



III. CALCULATION OF THE 
PHONON-ELECTRON HEAT TRANSFER RATE 
IN THE NORMAL AND SUPERCONDUCTING 
STATE 

In considering the electron-phonon thermal decou- 
pling, we assume that the electrons in the sample remain 
in local thermal equilibrium at a certain temperature T^i . 
We also assume that the phonons maintain local thermal 
equilibrium at a temperature Tph- We consider the possi- 
bility that the local electron temperature T^i differs from 
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the local phonon temperature, Tph and determine the re- 
sulting rate of heat transfer from the hot phonons to the 
cold electrons. This transfer rate is given bj*i^ 



dU 
Hi 



Kb -Lei 



(6) 



where the sum is taken over phonons of wavevector q, 
energy fuvqj^ velocity Cqj, and polarization j and n{x) is 
a Bose factor. All the electronic properties are contained 
within the sound attenuation aqj which is calculated sep- 
arately for the normal and superconducting states. 

The difference between the local temperatures of 
phonons and electrons is at most equal to the observed 
temperature difference between the two ends of the sam- 
ple (T' — T2), which corresponds to the entire range of 
temperature in the system. Since this temperature range 
is much smaller than the average temperature for these 
experiments, we can expand the above expression to first 
order in AT = Tph — Tgi, so that the heat transfer rate 
can be expressed as a electron-phonon thermal resistance 



(7) 



where T is the average temperature of the sample and 
AT can be thought of as the average temperature differ- 
ence between the phonons and electrons over the sample 
volume^. 

The fact that the thermal conductivity in PCCO is in- 
dependent of ficld^'' over the range from 8-13T suggests 
that the material is in the dirty limit, i.e. ujct << 1 
where is the scattering rate due to impurities and ujc 
is the cyclotron frequency. This suggestion is supported 
by the electrical conductivity data, which show a resid- 
ual resistivity corresponding to hT~^ /ks ~130K, which 
gives uJcT =0.1 for the highest field reported. The field 
has no effect on electronic transport in this limit, and we 
ignore it in our calculation of Rei~ph for this material. 
The condition for the low-temperature limit of the heat 
transfer rate: hr^^ >> {vf /cs)kBT (which corresponds 
to the low-frequency limit of sound attenuation for ther- 
mally excited phonons) is satisfiedi^i^iSii. With this, 
the well known result for the sound attenuation in met- 
als in the low- frequency limit is appropriate for PCCO. 
The normal state attenuation is given by 



2rjvno(fiwqj) 



< |/,(k,q)p (8) 



where is the impurity scattering rate obtained from 
the residual electrical resistance, uq is the density of 
states and M is the mass of a unit cell. The vector 
subscripts on the symbol, for the sound attenua- 
tion should remove possible confusion with the a used 
in Eq. |21to represent the electronic term in the thermal 
conductivity, so we will continue to use this standard 



notation. The electron-phonon interaction is described 
as in Ref. where is the coupling constant that 
sets the energy scale for the interaction and fj (k, q) is 
the dimensionless factor which describes directional and 
phonon- mode dependence. The explicit form of /^ fk, q) 
for each high-symmetry phonon mode is given in Ref. 1X81 
We determine the thermal resistance using Eqs. El d 
and|Sl Substituting Eq. [Slinto Eq. Eland using 



n[- — 

ksTph ' 

we obtain 



n{ -i^j 



knT 



keAT, (9) 



^ = kBATiksTf 



E 



< l/j(k,q; 



>FS 



qj 



dxx 



d^' 
(10) 



where p is the mass density and the flq integral is over all 
phonon directions (both /j(k, q) and Cqj depend on the 
mode and direction of the phonon but not on its energy) . 

The electron-phonon heat transfer rate in the normal 
state is thus seen to be proportional to T^ at temper- 
atures that satisfy hT~^ >> (vf /cs)kBT. This depen- 
dence can be understood by power counting as follows: 
two factors of temperature come from the Wq ^ depen- 
dence of the sound attenuation, three from the phonon 
density and an additional one from the energy of the 
phonon that appears in Eq. El with one power of tem- 
perature removed by the factor AT/T in Eq. El that 
describes the balance of energy exchange. We note that 
for the temperature T in Eqs. l^and lTUl we can use either 
Tph or Tel since these expressions are already first order 
in AT. 

By comparing Eq. [7|with Eq. ^| we obtain R^iLphi 
which is proportional to T^ , and use this along with the 
definition of Tb given in the Introduction to determine 
the downturn temperature. The result is 



. _ ll ^ 87rr(7)C(6)r^nofc|g^((|/,-(k,q)pc^5^) 



a \ 



{2TThfM 



(11) 

where the {{x)) indicates that the electron-phonon in- 
teraction has been averaged over the Fermi surface and 
phonon direction and summed over phonon modes. For 
the numerical estimates of Td described in Section IV 
we approximate the phonon averages by using the known 
expressions of /j(k, q) for high-symmetry directions and 
interpolating them to all other directions. The sound ve- 
locity of both transverse and longitudinal modes is taken 
to be isotropic. Transverse modes have a sound veloc- 
ity that is significantly smaller than that for longitudinal 
modes, which results in a heavier weight of the contribu- 
tion by the former to T^^. 

In the d-wave superconducting state, the low- 
temperature limit of Rei-ph, with its characteristic T^ 
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dependence is valid at temperatures such that 7 >> 
{vf /cs)kBT, where H~^j is the zero-frequency scattering 
rate of nodal quasiparticles by impurities. In this limit, 
the sound attenuation and hence R^i-ph are independent 
of 7, which is an example of the universal behaviour of 
transport properties in d-wave superconductorsiS^fi. The 
sound attenuation in this limit is related to the normal 
state expression of Eq. |S1 by 



qj 



l/.(kn,q)| = 



^ '"<|/,(k,q)P>p5 



(12) 



where k„ is the length of the wavevector to a node, and 
Hv2 is the slope of the gap at the node. We can use 
Eq. El in Eq. El to obtain the dependence of the 
electron-phonon heat transfer rate (note that the d-wave 
density of states is constant at temperatures ksT <^ j so 
that the power counting proceeds just as for the normal 
state). The value of T^'* for the superconducting state 
is obtained by multiplying the right side of Eq. 1111 by 
the factor a^Cqj/a^Cqj given by Eq. ^1 In the result- 
ing expression for Tjj'^, the matrix element |/j(k„, q)p is 
to be averaged over phonon direction and summed over 
phonon modes. 

Typical experimental values of 7 (which is related 

to Tat by 7 « 0.61y^Ao?iT^ if we assume unitary 

scattering)^ for most clean cuprate samples suggest that 
the strict low-T limit of Rei-ph occurs below the lowest 
temperatures used for thermal transport measurements. 
We thus must determine R^i^^f^, which becomes weakly 
dependent on 7, numerically over the experimental tem- 
perature range. The result of such a computation, using 
realistic values of 7 corresponding to between 1-1 5K as 
estimated from the position of the peak in k as a func- 
tion of temperature below Tc (see Ref. \Q and references 
therein), is well-described by a power law with an expo- 
nent n between 4 and 5. For example, by choosing 7 to 
correspond to a temperature of 3K, we find an exponent 
of n = 4.4, which can be used in Eq. When we com- 
pare our results to the in-field data on LSCO in Section 
IV, we will use the low-T (universal limit) result in order 
to simplify the comparison with the normal state result 
and the vortex state calculation described below. This 
result is only strictly valid within the rather narrow pa- 
rameter range of h(vf /cs)kBT « 7 << Aq, but should 
give a reasonable description for the case of LSCO over 
the temperatures of interest. If we were to allow 7 to 
vary over its experimental range, then we could adjust 
the value of n somewhat and perhaps improve the fit to 
the data. However, as seen in Section IV, the agreement 
between universal limit result and the available data is 
already very good. Thus in what follows we will use the 
universal limit result and will take fc^^7 = UK for LSCO 
(the zero-field result is independent of 7 but 7 does ap- 
pear in the field-dependent calculation below) in agree- 
ment with some estimates for this material from Ref. @. 
To study the magnetic field dependence of downturns 



observed in the vortex state of cuprates with a field ap- 
plied perpendicular to the Cu02 planes we follow the ap- 
proach of Volovik, and replace the frequency appearing 
in the superconducting Green's function with a Doppler- 
shifted value^^. The doppler shift is given by tw^ ■ k„ 
where Vs is the supercurrent velocity which depends 
on the position within the vortex-lattice unit cell, and 
k„ is a wavevector directed to the d-wave node. The 
nodal wavevector can be used instead of the true elec- 
tron momentum at the low temperatures of interest. 
The phonon-electron energy transfer rate can then be 
obtained as a position average over a simplified vortex- 
lattice unit cell^. Analytic results can be obtained in the 
clean, low-T limit defined by ■ k„ >> 7, {vf /cs)kBT, 
which can be expected to be valid in fields of several Tesla 
in many cuprates (here we refer to a typical value of -kn 
within the vortex-lattice unit cell). The formal similar- 
ity between the expression for the phonon attenuation 
rate and the thermal conductivity results in an identical 
field dependence for the electron-phonon heat transfer 
rate and the heat conductivity. The field-dependent at- 
tenuation in the superconducting state is given by 



AH) = (v, -k,, 



Tj^^V2k„ 



a; 



i(0) 



(13) 



where V2 is the magnitude of the slope of the supercon- 
ducting gap at the node and Q;^(0) is given by Eq. El 
The expression for the superconducting impurity scatter- 
ing rate in the clean, unitary limit 7 — t^^Ao/2vs • kn 
has been used. The electronic thermal conductivity is 
given in the same limit by 



T 



TAj^V2kn T 



«e;(0) 



(14) 



Ea.'sll3landll4lare to be ave rag ed over a unit cell of the 
vortex lattice. We follow Ref. 22 in performing the vortex 
unit cell average appropriate for fields applied perpendic- 
ular to the planes and assume unitary scattering. Both 
the thermal conductivity and the sound attenuation (and 
consequently R~i^_pfJ are thus found to be proportional 

to H/ Hc2- If we add the zero field universal result 
we have a decent description of the field dependences of 
a{H) and R^i^p^ for applied fields up to some fraction 
of Hc2- We then substitute into Eq. [SI to obtain 



Kel/T 



l+v{H)r(T/TD) 



with 



fl{H) = 1 + 



H a' 

Hc2 01 



(15) 



(16) 



where a' ~ a^Sir /8)\/ AQh~^TNa, and a is a factor of 
order unity that arises because we roughly approximated 
the vortex lattice geometry. The field dependence of the 
expression is contained entirely in ri{H). The factors r 
and Td are to be evaluated at = 0. 
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IV. COMPARISON WITH THE DATA 

In Fig. [21and Fig. Elwe show data for the field-induced 
normal state of optimally doped PCCOi and for super- 
conducting, optimally doped LSCO^, respectively. The 
plots are of the extracted electronic contribution to k/T 
versus T. Thus the phonon contribution has been sub- 
tracted from the data to leave only that of the electrons. 

According to our discussion in Section II, the phonon 
contribution could be identified by fitting the data at 
temperatures well above To to Eq. In practice, the 
second term in Eq. |S1 often does not adequately describe 
the phonon contribution (for example, slightly weaker 
temperature power-laws are often needed as discussed in 
Ref. 0). For our purposes, as long as the phonon con- 
tribution is determined from data at temperatures well 
above To, then once it has been subtracted we can use 
our results above and consider only the first term in Eq. 
El In Ref. Q the authors describe in detail how the 
phonon subtraction was performed for the PCCO data. 
For the LSCO data, we performed the phonon subtrac- 
tion by first fitting the data of Ref. at temperatures 
above the downturn (and below the maximum tempera- 
ture shown on Fig. 2)) to the form k/T = a-\- (3T™ where 
m was between 1 and 2 and then subtracting the (3T™ 
term. For both the LSCO data and the PCCO data, the 
electron contribution to k/T is seen to be independent of 
temperature for T » Tjy. 

We compare the result of Eq. in the normal state 
to the PCCO data in Fig. 13 For the result shown we 
used the parameters: a — 1.75mWK~^cm~-'^ that we take 
from the data at temperatures above the downturn, ?i = 5 
that we calculated in Section III, — 160mK, which 
was chosen to fit the data and r = 29 that was chosen 
to fit the data. Evidently Eq. O captures the low-T 
behaviour of the data. The large value of r indicates 
that heat is carried into the sample almost exclusively 
by phonons. Indeed, taking r = oo gives a reasonably 
good description of the PCCO data. This can be seen 
from Fig. O in which wc show the result of Eq. [S] for 
various values of r. 

In order to determine whether the value of To = 160 
mK is realistic for this material, we can compare this 
value with an estimate based on Eq. ^2 AH the parame- 
ters in Eq. II H are known for cuprates except the electron- 
phonon coupling constant, g^. The value of has 
been obtained from sound attenuation measurementsiS 
on Sr2Ru04, another material to which the given form of 
the electron-phonon interaction is applicable, and found 
as = (12eV)^. Although the electron-phonon cou- 
pling constant may vary significantly for different mate- 
rials, we can expect order-of-magnitude agreement be- 
tween the value of for Sr2Ru04 and that for PCCO. 
Using Eq. 1111 and the value To = 160 mK, we calculate 
the electron phonon coupling constant for PCCO and 
find g^ = (6eV)^. Thus, the coupling constant for PCCO 
determined using the fit above roughly agrees with the 
value measured for Sr2Ru04. 
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FIG. 3: Kei /T, normalized to its high-T value, determined by 
Eq. Ois plotted versus T for various values of the "contact" 
resistance r. 



In Fig. ^the result of Eq. Ofor the case of the d-wave 
superconductor in the low-temperature limit as a func- 
tion of applied magnetic field is compared to the LSCO 
data. Here, the zero-field result was obtained using the 
values: a = O.lSmWR-^cm-i, n ^ 5, To = 80mK, 
r = 0.86 that were chosen exactly as for the PCCO case 
described above. In order to obtain the 13T result, we 
use Eq. I15l with the values of a, n, Td, and r fixed by the 
zero-field fit (the two curves in Fig. ^were measured for 
the same sample). The remaining unknown factor of or- 
der unity, a, that appears in Eq. 1151 is obtained from the 
13T data at temperatures well above the downturn (we 
thus find a = 0.5). With this done, the low-temperature 
downturn for the 13T curve is determined without any 
adjustable parameters. The clear agreement with the 
data thus provides strong support for our model. 

To check whether the value of Td = 80 mK for su- 
perconducting LSCO is reasonable we compare it with 
the value determined in Section III. We note that the 
electronic thermal conductivity in the universal limit 
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FIG. 5: Main Panel: Kei/T, normalized to its high— T zero- 
field value, versus T in varying magnetic fields as given by Eq. 
1151 From the bottom, the curves are for fields of 0, 0.05, 0.1, 
and 0.2 Hc2 applied perpendicular to the Cu02 planes. Inset: 
The same quantity plotted versus H at various temperatures. 
From the bottom, the curves are for temperatures of 0.1, 1, 
and 10 Td. 



by any means other than via phonon-electron heating 
within the sample. We have not attempted to describe 
this complicated process and have treated r as an un- 
known parameter. However, if we assume that the main 
contribution to r comes from the resistance of the current 
contacts, which is usually measured during thermal con- 
ductivity experiments, then we can compare the values 
of r obtained above to the measured values. We convert 
the thermal conductivity of the sample measured well 
above the downturn to an electrical resistance by using 
the Wiedemann Franz law, according to which 



Re 



(17) 



where A is the cross-sectional area of the sample and Lq 
is the Lorenz number. For the PCCO measurements, the 
contacts were reported^ to have a measured resistance of 
roughly 1 fl. So, using the reported sample dimensions. 
Is — 1 mm, A = S8 fim x720/im and using Rei(c) = 
(the factor of two mentioned in the discussion following 
Fig. ^has been included) we get r — 240. This is sub- 
stantially larger than the best- fit value of 29, but from 
Fig. Ois seen to still be reasonably consistent with the 
data. For the LSCO data, the reported contact resistance 
was roughly 10 mfl. The dimensions of the sample24 for 
which the data above was obtained were Is = 0.97 mm, 
A — 1.26mmx0.212 mm. From this we obtain r = 1.8, in 
rough agreement with the value r = 0.9 used above. We 
do not expect precise agreement between the r obtained 
above and the measured contact resistance, but it is clear 
that the values of r used in our fits are not out of line 
with the measured values. 

The data shown above was chosen for being represen- 
tative of that in which downturns are seen. Other ex- 
amples of downturns, some of which were discussed in 
Section I, typically resemble those shown above and oc- 
cur at roughly the same temperature. 



DISCUSSION 



is related to the normal state Boltzmann expression 
by an expression similar to Eq. 1121 namely Kq /T — 
/T{t^^ /Trv2kn)- As a result, the value of Td de- 
termined using the low-T results for both Rei-ph and 
a = kq/T is roughly the same in the superconducting 
state of cuprates as it was for the normal state. From 
Eq. ^](also, recalling the discussion immediately follow- 
ing Eq. 1121) we find that the value of Td = 80mK corre- 
sponds to — (7eV)^. The agreement of the electron- 
phonon coupling constants for PCCO and LSCO is ex- 
pected since is a normal state property of the Cu02 
planes. 

The parameter r = Rei{c) / Rei determines the extent 
of the low-T downturn, i.e. the value of k/T aX T — Q. 
The value of r corresponds to the thermal resistance to 
heat being transmitted to the electrons in the sample 



We have shown that the low-T downturn in the elec- 
tronic contribution to thermal conductivity, for PCCO 
in the normal state and LSCO in the superconducting 
state, results from the loss of thermal equilibrium be- 
tween phonons and electrons caused by poor heat transfer 
between these systems at low-temperature. The electron- 
phonon heat transfer rate has been shown by calculation 
to have both a characteristic temperature dependence in 
the normal state, and a characteristic temperature de- 
pendence and magnetic field dependence in the super- 
conducting state. The fact that these dependences are 
in excellent agreement with experimental observations in 
PCCO and LSCO and provides a striking confirmation 
that poor electron-phonon heat transfer at low tempera- 
tures is indeed the right mechanism to explain the down- 
turns in these materials. Clearly this is a mechanism 
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that should also be considered in connection with ther- 
mal conductivity downturns in other materiak^^^®. If 
one seeks to measure the electronic contribution to k/T 
at zero temperature, the role of electron-phonon thermal 
decoupling has to be carefully considered. Alternatively, 
one could study the downturns to extract information 
about the electron-phonon interactions, without regard 
to the underlying thermal conductivity. 

Thermal conductivity data at temperatures below To 
that shows a significant downturn could be used to study 
electron-phonon interaction effects in unconventional su- 
perconductors and other systems of interest. In d-wave 
superconductors, the electron-phonon heat transfer rate 
calculated within the dirty d-wave approach has many 
of the same features as the thermal conductivity. We 
showed that it has the same magnetic field dependence 
at low-temperature, and also shows universal behaviour 
(albeit at a temperature which is smaller by a factor 
of the ratio of the speed of sound to the fermi veloc- 
ity, Cg/vf than that at which universality occurs in the 
thermal conductivity>i22S2i) At higher temperatures, the 
field and temperature dependences of the heat transfer 
rate have different forms from those of thermal conduc- 



tivity, and provide complementary information about the 
electrons. The basic form of the electron-phonon matrix 
elements, which affects the heat transfer rate, also has 
unusual properties in layered square-lattice tight-binding 
materialsiSS^. 
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